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Thermoelectric energy harvesters can have a much higher conversion effi¬ 
ciency by implementing quantum dots/wells between the high temperature 
region and the low temperature region. However they still suffer a limitation 
of the maximum output power, represented by the maximum Ai? (maximum 
energy gap of two quantum dots/wells layers). In this work, we use the 
piezoelectric material in the high temperature region, which has conceptually 
addressed the problem of the maximum power limitation. Full analysis of 
device physics including comparison with the existing technology and quantum 
simulation has been conducted to validate this concept. Results show that with 
the new concept, the maximum output power has been increased by at least an 
order of magnitude with the same power input and identical device dimensions. 
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I. INTRODUCTION 


Thermoelectric energy harvesting has the re- 

Q rse effect as the thermoelectric refrigerator 
, and has rbe^ ^udied extensively in re¬ 
cent decades (21 [Sj , alongside oth^ 110 ,- 

^ble nanometre sized energy harvesters |6| |7| 
[sf. With regards to the recent development in 
thermoelectric harvesting, several literatures re¬ 
ported that implementing quantum dots or wells 
into the dpcice-can-^ignihcantly increase device 
efficiency (sl 1^ llj, attributed to the energy 
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hltering effect of the quantum dots/wells, whi4S|h 
coincides with what has been suggested in [51. 
Theoretical investigation of quantum dots based 
thermoelectric harvesters ha^been conducted, 
which was reported in Ref. |9| that the maxi¬ 
mum scaled output power with other parameters 
being optimized appears at around AE=6kBT, 
AE being the difference of the energy levels of 
two quantum dots layers on the left and right 
sides of the central cavity. Question is then 
whether there is any solution to overcome this 
limit. Right from the invention of the seebeck 
and Peltier effects up to now, people have been 
using doped semiconductor materials as the cen- 
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tral cavity with the aim of increasing the electri¬ 
cal condnctivity and reducing thermal conduc¬ 
tivity for a higher hgure of merit. Specihcally 
for the quantum dots thermoelectric harvesting, 
very large energy level difference (larger than 
the optimized AE) of quantum dots in left and 
right sides will cause electrons difficult to en¬ 
ter/exit the central cavity (schematic illustration 
in Fig. d]). It is seen from the schematic diagram, 
the chemical potential of the central cavity gen¬ 
erally deployed has been a flat line. Changing 
the flat potential of the central cavity to a tilted 
shape will be a conceptual advancement, as ide¬ 
ally the left end of the chemical potential of the 
central cavity should be align with the energy 
level of the left quantum dot, and the right end 
should be align with the right side quantum dot 
in order to eliminate reflections at the staircase 
in the potential diagram. It is worth noted that 
the new concept applies to increasing the effi¬ 
ciency of the thermoelectric cooling as well. Out 
of many semiconductor materials, a poled piezo¬ 
electric material has a ramped chemical poten¬ 
tial, which perfectly meets this application. In 
this article, a piezoelectric material is used for 
the central hot region. Theoretical investigation 
conducted here shows that a signihcant increase 
in energy efficiency has been observed. Detailed 
analysis is described and presented in below sec¬ 
tions. 


II. QUANTUM DOTS 
THERMOELECTRIC HARVESTER 

In order to understand the motivation of this 
concept, the optimization of the previously re¬ 
ported quantum dots thermoelectric harvester 
is demonstrated in this section. As schemat¬ 
ically shown in Figs. [I](a) and 1(c), electrons 
entering to the central hot area from left side 
electrons bank gain energy from phonon-electron 
interaction, and energized electrons flow to the 


right side of the electrons bank forming a elec¬ 
trical current loop, which operates as a battery. 
Fig. [IKc) schematically displays that electrons 
having lower energies than the central potential 
will be precluded from entering the central cav¬ 
ity, which is thought the underlying reason of 
the device being limited by the maxirp|im energy 
gap AE. Following the procedure in |9|, the out¬ 
put power of the harvester can be derived from 
the conservation laws for electrical charge and 
energy 


II - fcL + Ir- fcR = 0 

T, - ELifL - Icl) + - fcR) = 0 ( 1 ) 

27 

The hrst part of the Eq. ([ 1 ]) represents that the 
total electrical charge equals to zero, which in¬ 
dicates that the charges flowing from the left 
electrode reservoir to the central cavity equals 
to the charges flowing from the central cavity 
to the left side reservoir. Where /l, fci, Jr, fcR 
are Fermi Dirac distributions characterizing the 
cavity’s occupation function, which are fcL = 
f{EL - ficTc), h = f{EL - fcR = 

/{Er - /ic,Tc), fR = /{Er - fiR,TR). The 
above expressions represent the occupation den¬ 
sity of electrons within the energy band (band¬ 
width is 7 ) centred about each energy level. The 
general form of the Fermi Dirac distribution is 
f{E - ^,T) = 1/(1 + second 

part of the Eq. ([H) describes that the total energy 
(external heat energy, heat current from the left 
to the right reservoirs) is zero, h is the Planck’s 
constant. Re-organize the Eq. CDi, one can get 
the heat current J 


J = 


2-/AE 

h 


ifcR - Jr) 


( 2 ) 


The electrical current / is 

J _ _e^ _ ^ejjfcR - /fl) 
~ AE ~ h 


( 3 ) 
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Normal material 



(a) 




FIG. 1. (a) Schematic illustration of the previously reported energy harvester based on quantum dots. Central region is 

connected to a heat source, where energy of electron has been increased through phonon-electron interaction, (b) Proposed 
new piezoelectric thermoelectric harvester. Poled piezoelectric material has been used in the central region, so that a chemical 
potential slope fj,p is present. The chemical potential slope of the piezoelectric layer has signihcant impact on the maximum 
output power analysed in the text, (c) 3D illustration of (a) in which electron having lower energy than the central energy level 
is bounced back, (d) Diagram of the new concept where much more electrons with lower energies can proceed to the central 
cavity, leading to much improved energy output. 


where e is the charge of an electron, and AE = 
Er — El- The efficiency of the system de¬ 
fined as the ratio of the harvested power P = 
\{^L — AiR)-^|/e to the heat current J. are 

chemical potentials of the left and right elec¬ 
trodes. The efficiency rj is then expressed as 
V = (l/^T “ It is then derived that 

the /i = l/iL — fiR\ < AE, as the 17 < 1. The 
output power is expressed as 


P = -^{fcR - fn) 


(4) 


results are shown in Fig. [2l For the parame¬ 
ters set in the figure caption, the maximum P 
is at around 6k rT, at which the optimized /i is 
around 0.43Ai?, and the 7 = ksT. The above 

a suits match closely with what was reported in 
]. With the question of increasing the maxi¬ 
mum power, it is envisaged that a ramped chem¬ 
ical potential may result in a better performance, 
as electrons with lower energy can have the op¬ 
portunity to enter the central region (Fig. [T]^d)). 
Device physics of the new configuration is de¬ 
scribed in the following section. 


Numerical simulation can be performed to arrive 
at the maximum output power as the function of 
the AE. In the simulation, and 7 have been 
optimised to have the maximum P. Simulated 
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FIG. 2. (a) Pma^h/iksTf versus AE/ikaT) for AT = T 

and is optimized, and the level width 7 is chosen as kaT. 
AT = Tc — Tr and T — {Tc + Tr)/2. (mr^l = ±/i/2 + {El + 
Er)I 2, He = {El + Er)I2. hc is the chemical potential of 
the central region, (b) AE is fixed at the optimized value 
5.6AE, scaled maximum power Pmaxh/{kaT)^ as functions 
of H and 7 . 


III. DEVICE CONCEPT 
INCORPORATING PIEZOELECTRIC 
MATERIAL 

In the new device architecture (Figs. [11(b) and 
1 (d)), with the aim of overcoming the limitation 
of the maximum output power obtained from the 
previous analysis (~6 /cbT), a piezoelectric mate¬ 
rial has been chosen to make the central cavity. 
For a poled piezoelectric material, the chemi¬ 
cal potential is no longer a flat line, instead a 
slope [12|, whose gradient is dependent on the 
material properties. The potential gradient of 
the poled piezoelectric materials is understood 
from aligned electrical dipoles. A piezoelectric 


material can have both electrical conductivity 
and piezoelectricity simultaneously, for example 
bismuth titanate (BIT) [isj. Moreover due to 
the high transition temperature, BIT ceramics 
are good candidates feauhigh temperature piezo¬ 
electric applications [1^. Well processed piezo¬ 
electric material will meet the requirement for 
the hot cavity in the thermoelectric harvester in 
the viewpoint of practical realization. The en¬ 
ergy levels of the quantum dots can be tuned 
by the external bias, alternatively they can be 
re-conhguredJay the externally applied mechan¬ 
ical stress 1^ Ig. In the new device, the 


chemical potential of the central region becomes 
= {El + Eji)/2 ± /ip/ 2 , where the are 
the upper and lower ends of the chemical poten¬ 
tial of the piezoelectric layer, /ip is the difference 
between the upper and lower ends of the chemi¬ 
cal potential of the piezoelectric material. Hence 
the fcR becomes 


exp[(A£; - Hp)/{3kBT)] 

Substituting the Eq. [5] to the Eq. 01 numerical 
analysis has been conducted using optimized n 
and 7 . The /ip independent of AE has been 
initially designated to hxed values as the func¬ 
tion of fcpT. The results (shown in Fig. [3]) dis¬ 
play a signiheant increase of the scaled maximum 
power Pmaxh/{ksTY from around 0.4 to 5 with 
the /ip increasing from 0 to SksT. However it 
does not demonstrate the similar scale of the in¬ 
crease for the optimized AE, which can be ob¬ 
served to increase from around b.GksT to around 
12/cpT as /ip increasing from 0 to 15/cpT. This 
is due to that the /ip is given hxed values while 
the AE varies. Further investigation is there¬ 
fore performed with /ip being dependent upon 
the AE. In the analysis, the /ip has been des¬ 
ignated as a ratio of AE, i.e. from 1/5 of AE 
to 4/5 of the AE. Numerical calculation (shown 
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FIG. 4. (a) Maximum output power versus energy band gap 

of quantum dots for various piezopotentials proportional to 
AE. (b) Optimized AE versus the ratio of fip to AE. 


FIG. 3. (a) Calculated maximum output power 

Pmaxh/{ksT)^ relating to energy gap between two layers of 
quantum dots AE/iksT) for various fixed piezopotentials 
fip/iksT). (b) Optimized AE in relation to the piezopo¬ 
tentials fj.p. 

in Fig. Sj) demonstrates that apart from the ob¬ 
vious increase of the maximum power, the opti¬ 
mized AE has also increased to around 80k bT 
for /tp = 0.9Ai?. Extended analysis has been 
conducted, which shows that when the pp ap¬ 
proaches to AE, there is no maximum AE, in¬ 
stead the maximum power continues increasing 
with the increasing AE. 

IV. QUANTUM ANALYSIS 

Further analysis is needed using quantum me¬ 
chanics theory to reveal the physics behind the 


performances of the new proposed device. As 
the fundamental element of the device is a qi^.- 
tum dot, quantum theory has to be used [17 1. 
The core of the quantum theory is based on 
the Schrodinger equation, here both the time- 
dependent and time-independent Schrodinger 
equations are used in the simulation, which are 
shown in below 

+ V{x,t)]ilj{x,t) (6) 

where 'ip{x, t) is the wave function of an electron 
as the function of time t and location x, h is 
Planck’s constant divided by 27r. m is the mass 
of the electron and V{x, t) is the potential di¬ 
agram through which the electron passes. The 
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FIG. 5. Time-independent analysis. To simplify the calcu¬ 
lation, let h=l, m=l, length in x direction is 100. Ai5=14, 
fj,p — AE, p=10. (a) Ip versus position in x and the energy 
diagram of device with a standard material, (b) ip versus posi¬ 
tion in X and the energy diagram of device with a piezoelectric 
material. 


time-independent Schrodinger equation simpli¬ 
fied from the Eq. |6] is 

EiP{x) = + y{x)]ip{x) (7) 

where only the static state of the probability 
along the location x of an electron with certain 
energy level B can be calculated. One dimen¬ 
sional numerical simulation has been conducted 
using Eqs. [6] and [71 Results from calculations 
have been shown in Figs. [5] and El Fig. El shows 
the calculated distribution of the wave function 
ijj for the electron with the initial energy level 
of 8.8 a.u. along both the energy diagrams rep¬ 
resenting the device with a traditional material 


and the device with a piezoelectric material re¬ 
spectively. All parameters used in the simulation 
are defined in the figure caption of the Fig. El It 
can be summarized that for the standard mate¬ 
rial there is very little wave penetrating to the 
central region as the energy level of the central 
material is 10, which is higher than the elec¬ 
tron energy level. However for the device con¬ 
structed by a piezoelectric material, there is a 
good portion of wave passing to the central re¬ 
gion as highlighted by the red elliptic shape in 
Fig. Etb). It should be noted that the amount 
of wave through to the central region will be en¬ 
ergized by the external heating source and even¬ 
tually going to the left side of the device form¬ 
ing an electrical current loop with external loads. 
While for the device with the normal material, 
electrons with less energy than the central en¬ 
ergy level will die out immediately. Fig. El shows 
the dynamic simulation as to an electron pass¬ 
ing through the two potential diagrams using 
the one dimensional time-dependent Schrodinger 
equation. At a sequence of time points, the 
wave packet travels along the axis x for the two 
different energy diagrams have been calculated 
using simplified parameters. At the beginning, 
two cases are pretty much similar. A little por¬ 
tion of the wave packet has been reflected at 
the potential step of the device with the cen¬ 
tral cavity made of a traditional material. In 
comparison, it is a very smooth transition with¬ 
out any reflection on the potential diagram of 
the piezoelectric material. As the time develops 
further it can be seen that in the normal mate¬ 
rial, the wave packet travels faster than in the 
piezoelectric material, leading to less phonon- 
electron interaction time. It is seen from the 
results that three advantages are demonstrated 
by the novel piezoelectric quantum dots thermo¬ 
electric energy harvester: 1), more electrons with 
lower energy can reside in the central region, 
evidenced by the time-independent Schrodinger 
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FIG. 6. Time-dependent quantum analysis. To simplify the calculation, let h=l, m=l, length in x direction is 100. AI5=14, 
= Ai?, /i=10. (a)-(e) moving wave packet from the left to the right, energy level of electron is 8.8 a.u., at time points of 0, 
2000, 4000, 5000, and 6000 respectively. This is for the device with a piezoelectric material, (f)-(j) moving wave packet of the 
electron (E=8.8 a.u.) for the device with normal material in the central. Starting position of both cases is from x=10, initial 
amplitude of both the wave packets is set to 1, initial width of the wave packet is 20, Aa;=0.1, and At=0.003. 


analysis. 2), less reflection can be achieved, sup¬ 
ported by the time-dependent Schrodinger anal¬ 
ysis. 3), slower wave packet travelling speed of¬ 
fering longer thermoelectric transfepring time is 
observed. It was described in the |9| that with 
the accumulative effect taken into account, the 
quantum dots thermoelectric energy harvester 
could reach the value of O.IW per cm^ at temper¬ 
ature difference of IK. From above analysis , the 


new piezoelectric material based quantum dots 
thermoelectric harvester could have the maxi¬ 
mum power output of more than IW per cm^, 
which presents the highest efficiency for the en¬ 
ergy harvesters of this kind. 

V. CONCLUSION 

To conclude, with the objective of achieving 



























































higher power efficiency of the quantum dots ther¬ 
moelectric energy harvester, a piezoelectric ma¬ 
terial is used to replace the standard material 
in the high temperature region. According the 
analysis, the power output of the new device will 


be at least 10 times higher. Quantum mechan¬ 
ical analysis reveals that with a tilted chemical 
potential of the piezoelectric layer, the electron 
passes through the central region with minimal 
reflection, and more electrons are energized by 
the thermal energy in the central hot region. 
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